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Abstract 

By applying the recently developed Loop Regularization(LR) with string-mode regulators to supersym- 
■ metric field theories, we explicitly verify the supersymmetric Ward identities in several supersymmetric 

models at one-loop level. It is interesting to observe that supersymmetry is a so remarkable symmetry 
that the supersymmetric Ward identities hold as long as a regularization scheme is realized in the exact 
four dimensional space-time with translational invariance for the momentum integration, and the gauge 
symmetry can be maintained once the regularization scheme preserves supersymmetry and satisfies the 
consistency condition for logarithmic divergences. As a manifest demonstration, we carry out a complete 

CI 

one-loop renormalization for the massive Wess-Zumino model by adopting the LR method, it is found that 
all the quadratic divergences cancel out and the relations among masses and coupling constants hold after 
renormalization, which agrees with the well-known non-renormalization theorem. It is concluded that the 
Oh| LR method preserves not only gauge symmetry but also supersymmetry. A simple and definite derivation 

O ■ of Majorana Feynman rules is found to be very useful. 
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I. INTRODUCTION 



Supersymmetry has attracted physicists for several decades since it was proposed in 1970s[l|]. 
As it is well-known that symmetry has played an important role in particle physics, and three of the 
four basic forces in nature are governed by gauge symmetries and have successfully been described 
by quantum field theory. While quantum field theories are bothered by the infinities which must be 
regularized to be well-defined. On the other hand, whether the symmetries of classical Lagrangian 
still hold in the quantum level remains an important issue, this is because some times it is difficult 
to distinguish between a real anomaly and an apparent violation of the symmetries due to the use of 
a symmetry-violating regularization method. In general, when a symmetry of original Lagrangian 
is still a symmetry of a full quantum effective action, such a symmetry is regarded to be preserved 
in the quantum level, but there are several exceptions such as chiral anomaly. Thus one may ask 
whether supersymmetry is a symmetry of the full quantum theory. This question has been studied 
in a regularization- independent way in ref.Q], and the answer is yes. This means when investigating 
the quantum effects of the supersymmetric theories, one must adopt a supersymmetry-preserving 
regularization method. 

Several regularization methods have been applied to supersymmetric theories, such as dimen- 
sional reduction(DRED)[3], differential regularization [3] and the so-called implicit regularization 0], 
among them DRED is the most common one. It has been shown that DRED can preserve super- 
symmetry in several models^, Q, Sj- Strictly speaking, DRED is mathematically inconsistent (9I. flo| 
to be applied to the supersymmetric theories, which is similar to the case when it is applied to 
the chiral theories, this is because both supersymmetry and the definition of 75 require an exact 
dimension. A consistent regularization method that can be applied to all possible cases in quantum 
field theories is needed. In this sense, the recently developed loop regularization(LR) with string- 
mode regulators 11, may deserve a special attention, it has successfully been applied to the 
calculations of triangle anomaly of QED with clarifying the possible ambiguities caused by 75 13], 
the evaluation of a consistent coefficient of the CPT and Lorentz symmetry breaking Chern-Simons 
termjlil]. the computation of all the one-loop renormalization constants for the non-Abelian gauge 
theory and the determination for the coefficient of QCD (5 function [l5|], and the derivation of 
the chiral effective field theory with a dynamically generated spontaneous symmetry breaking [3]. 
The key concept of this new regularization method is the introduction of the irreducible loop 
integrals (ILIs) which are evaluated from Feynman integrals by using Feynman parameter method. 

It has been shown that the LR method can preserve the non-Abelian gauge symmetry, and 
meanwhile maintain the divergent behavior of original field theories. In particular, the LR method 
is realized in the original four dimensional space-time with translational and Lorentz invariance 
even if two intrinsic mass scales are introduced, thus it can balance the bosonic and fermionic 
degrees automatically and there is also no ambiguity about the definition of 75. It is then believed 
that this method will preserve supersymmetry as well. In this paper, we will investigate the 
applicability of LR method in supersymmetric theories. 

The paper is organized as follows: in section II, we briefly introduce the symmetry-preserving 
loop regularization with string-mode regulators. In sections IlIII and ITVl we will verify the supersym- 
metric Ward identities for the massless Wess-Zumino model and massive Wess-Zumino model 17] 
separately and show that the LR method indeed respects the Ward identities. As Ward iden- 
tity is the reflection of symmetry in quantum level, we then arrive at the conclusion that the LR 
method is also a supersymmetry-preserving regularization for supersymmetric models. In section 
PVT we consider the super Yang-Mills theory as the testing ground to explicitly demonstrate the 
supersymmetric Ward identity and show that the LR method does preserve supersymmetry, and 
meanwhile the gauge symmetry is maintained only requiring the consistency condition for logarith- 
mic divergences. Note that the conventional dimensional regularization was shown to break the 
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Ward identity in such a model[6j, thus an alternative check by using LR method in our present 
paper is nontrivial. In particular, we will demonstrate that as long as the Dirac algebra for 7 
matrices are carried out in four dimensional space-time, and the shift of integration variable can be 
safely made, the supersymmetric Ward identities are preserved, which is actually independent of 
any concrete prescription of regularization methods. Namely, as long as the regularization scheme 
is realized in four dimensional space-time with translational invariance for momentum integrals, 
like the LR method, it then preserves supersymmetry. In section PVTl as an explicit demonstration, 
we will carry out the one-loop renormalization for the massive Wess-Zumino model by using the 
LR method , and all the obtained results agree with the well-known non-renormalization theorem. 
Our conclusions and remarks are given in the last section. The detailed derivation of Majorana 
Feynman rules is presented in the appendix. 



II. SYMMETRY-PRESERVING LOOP REGULARIZATION 

It has been shown in 11, that all one loop Feynman integrals can be evaluated into the 
following 1-fold ILIs by using the Feynman parameterization method: 
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d 4 k 1 
(2^) 4 (k 2 -M 2 ) 2 + a ' 



I-2a»u = J (2 7T) 4 (fc 2 -M2)3+a' " = -1,0,1,2,... 

r d 4 k kp~k v kpk a 

l-2a »vpa ~ J ( 27r )4 ( fc 2 _ M 2)4+a V' L ) 

with I2 and Iq corresponding to the quadratic and logarithmic divergent integrals. Where the 
effective mass factor M 2 is a function of the external momenta pi , the masses of particles m, and 
the Feynman parameters. 

In general, the loop momentum independent M 2 can be extended to include linear term in k, 
which can be understood as a part of the definition of the ILIs in the LR. The reason is as follows: 
Let M 2 {k) has the following general form including linear term in k 

M 2 (k) = M 2 + 2xk.p 

with x an arbitrary parameter. Then 

k 2 + M 2 (k) =k 2 + 2xp.k + M 2 = (k + xpf + M 2 - p 2 = (k + xpf + M 2 {p) = k' 2 + M 2 (p) 

with M 2 (p) = M 2 -p 2 which becomes independent of k, and k! = k-\-xp via translational invariance. 
Again the only thing must be paid attention is that one must follow the definition of ILIs to cancel 
out the k 2 in the numerator before regularization. 

When the regularized 1-fold ILIs satisfy the following consistency conditions Hi, LL2J : 



^2pu — 2 5 7^ I2 ' ^2[ivpa — gidfiudpa + 9^p9vu + 9^9 pv) 1 2 1 

lopv = l9nv Iq 1 } loiivpa = ~^\9pv9po + 9pp9va + 9po9pv) Iq- (2-2) 

the resulting loop corrections are gauge invariant. Here the superscript "R" denotes the regularized 
ILIs. 

Note that the introduction on the concept of irreducible loop integrals (ILIs) is crucial in the loop 
regularization 11, [13], where it has been shown that all Feynman loop integrals can be evaluated 
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to be expressed by the ILIs. From the definition of ILIs, one of the important properties is that 
there should be no k 2 in the numerator of loop integration, all the ILIs can be classified into the 
scalar type ILIs with the following loop integration 

1 



(k 2 - M 2 Y 

and the tensor type ILIs with the following loop integration 

k^k v • • • kp 
(k 2 - M 2 ) a 

In evaluating the Feynman loop integrals into ILIs, one should always perform the Dirac algebra 
and Lorentz index-contraction firstly to obtain the ILIs defined by the above "simplest" forms for 
the one loop case (for two loop and higher loop case, see ref.[ll]). Therefore, for the integration 



which should not be written as 



but it must be expressed as 



<T ■ k,k v /{k 2 - M 2 ) 2 
<T • hp* 



k 2 /{k 2 - M 2 ) 2 

then rewriting the k 2 in the numerator into (k 2 — M 2 ) + M 2 so as to cancel out the first term by 
the denominator. Thus the above Feynman loop integration is regarded to be evaluated into the 
ILIs and is given by the following form before regularization 

sT ■ kpKlik 2 - M 2 ) 2 = I 2 + M 2 * Jo 

From the above illustration, it is seen that in the spirit of "irreducible loop integrals" (ILIs), 
the integration 

gT ■ k,k v /{k 2 - M 2 ) 2 

is not an ILI, one should not regularize such a loop integration in the loop regularization method. 

A simple regularization prescription for the ILIs was realized to yield the above consistency 
conditions, its procedure is that: Rotating to the four dimensional Euclidean space of momentum, 
replacing the loop integrating variable k 2 and the loop integrating measure / d A k in the ILIs by 
the corresponding regularized ones [k 2 ]i and J[d A k]i: 

k 2 -► [k 2 ]i = k 2 + M 2 , 

N 

[d 4 k^ [ [d%t = lim V cf [ d 4 k (2.3) 

J J N,Mf ; _ J 

where M 2 (I = 0, 1, • • •) may be regarded as the regulator masses for the ILIs. The regularized 
ILIs in the Euclidean space-time are then given by: 

' d^k 



I R 2* = <(-!) , 



Y \ ^ n f ® k 1 

tv!m^ C/ J (2vr) 4 (fc 2 + M 2 + M2)2+«' 

^ - lim g f J ^g^ , « = -1,0, 1,2, ... 
tR - A(_-\\a \^ r N f d4k kpkykpkg 

i-Ta^pa ~ H 1J S^fl J (2vY {k 2 + M 2 + M 2 Y+« { > 



i 1=0 
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where the coefficients cj^ are chosen to satisfy the following conditions: 



N 

Jfc\/r2\n 



\im}2 c ?( M n n = (» = 0, !,■••) (2-5) 



with the notation lim^y^ denoting the limit lim N m^—hx' One may take the initial conditions 
-^o = Ms = and Cq' = 1 to recover the original integrals in the limit Mf — > oo (/ = 1,2, • • • 
) . Such a new regularization is called as Loop Regularization (LR) ll|, [12J . The prescription in 
LR method is very similar to Pauli-Villars prescription, but two concepts are totally different as 
the prescription in the loop regularization is acting on the ILIs rather than on the propagators 
in Pauli-Villars scheme. This is why the Pauli-Villars regularization violates non-Abelian gauge 
symmetry, while LR method can preserve non-Abelian gauge symmetry. 
As the simplest solution of eq. (|2.5|) . taking the string- mode regulators 

Mf = £ + IM 2 R (2.6) 

with I = 1, 2, • • •, the coefficients are completely determined 



Here Mr may be regarded as a basic mass scale of loop regulator . It has been shown in [Tj 
that the above regularization prescription can be understood in terms of Schwinger proper time 
formulation with an appropriate regulating distribution function. 

With the string-mode regulators for Mf and cf in above equations, the regularized ILIs 1^ and 
Iq can be evaluated to the following explicit forms 



11, 121: 



—i M 2 a 2 

® = l6^ {Mc2 -^ [/n #- 7 - + 1 + y2( W )]} 

T R i r, M c , M 2 m 

J ° = 16^ [/n /-^ + yo( Mf )] 



with (i 2 = (i 2 + M 2 , and 

N N 

lw = lim{ ^cf lnZ + H^cf /In/ }} = lE = 0.5772---, 

1=1 i=i 

y (x) = / da , yi{x) = 

Jo a x 

2/2O) = Vo{x) - Vxix), limyi(x) 0, i = 0,1,2 (2.9) 

x— >0 

AT 

M 2 = lim Mn V cFil In Z) = lim Ml I In V 

which indicates that the (i s sets an IR 'cutoff' at M 2 = and M c provides an UV 'cutoff'. For 
renormalizable quantum field theories, M c can be taken to be infinity (M c — > 00). In a theory 
without infrared divergence, (i s can safely run to (j, s = 0. Actually, in the case that M c —* 00 
and (a s = 0, one recovers the initial integral. Also once Mr and ./V are taken to be infinity, the 
regularized theory becomes independent of the regularization prescription. Note that to evaluate 
the ILIs, the algebraic computing for multi 7 matrices involving loop momentum }f, such as ^7^^ 
should be carried out to be expressed in terms of the independent components: 7^, a^ u , 75 7 M , 75. 
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We shall directly show that loop regularization is manifestly translational invariant in spite of 
the existence of two energy scales, which is a very important feature in applying to super symmetric 
theories in this paper. To see that, we shall verify that the regularized ILIs should arrive at the same 
results whether the loop regularization prescription is applied before or after shifting the integration 
variables for momentum. For an explicit illustration, let us examine a simple logarithmic divergent 
Feynman integral: 



L = I 7^113 H2 7I ^2 ^2 (2-10) 



d 4 k 1 1 
(27r) 4 k 2 — m 2 (k — p) 2 — m\ 

As the first step of loop regularization, we shall apply the general Feynman parameter formula 



1 r(a!i + • • • + ol ) f Xl r%n-2 

ai an 57 = ~F? A TV T" dx l dx 2 ■ ■ ■ dx n -i 

a^a^-'-an" T(qi) • • • T(a n ) Jo Jo Jo 

(1 - xi) a ^ l { Xl - x 2 ) a ^ 1 - • -x^I 1 
[ai(l - x x ) + a 2 (xi - x 2 ) + • • • + a n x n _i] a i+-+«* 

to the Feyman integral and obtain the following integral 

L=i£ir dx 



(2.11) 



(2vr) 4 Jo {(1 - x)(k 2 - m\) + x[(k - p) 2 - m 2 }} 2 
d 4 k f 1 , 1 



dx- 



(27r) 4 jo {(k — xp) 2 — [(1 — x)m\ + xm 2 — x(l — x)p 2 ]} 2 
f 1 f d 4 k 1 
= io dX J (2tt) 4 [(k - xp) 2 - M 2 ] 2 (2,12) 

with M 2 = (1 — x)m 2 + xm 2 — x(l — x)p 2 . 

By making Wick rotation and applying the loop regularization prescription before shifting the 
integration variable, i.e., rewriting the momentum factor {k — xp) 2 into (k — xp) 2 = k 2 — 2xp.k+x 2 p 2 , 
then replacing k 2 by k 2 + M 2 , namely 

(jfe - xp) 2 = k 2 - 2xp.k + x 2 p 2 -> k 2 + M? - 2xp.k + x 2 p 2 = (k - xp) 2 + M? (2.13) 

we then obtain the regularized Feynman integral 

LR = 1 g ^ I ^ / Wf [(k-xp) 2 + M 2 + M 2 ] 2 (2 ' 14) 
which becomes a well defined integral, so that we can safely shift the integration variable: 

L R = dx lim Vcf / /n , dn9 -= r-joTn (2.15) 
io ^Mf^ 1 7 (27r) 4 (fc 2 + M 2 + M 2 ) 2 y ' 

The same result can be arrived by using the standard procedure of loop regularization with first 
shifting the integration variable for momentum, which yields the standard scalar type ILI 

rl r d A k 1 I' 1 

L o = / dx 7^1779 = / dxI (2-16) 



(2vr) 4 (k 2 - M 2 ) 2 Jo 
after applying the loop regularization prescription, the same form is reached 



L "^[ d '^I^WTjhMff sLR (2 ' 17) 
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which shown that in loop regularization method, one can safely shift the integration variables and 
express all the Feynman integrals in terms of ILIs before applying for the regularization prescription. 

From the above explicit demonstration, it is seen that the loop regularization is indeed trans- 
lational invariant. In fact, this property also allows us to eliminate the ambiguities and make a 
consistent calculation for the chiral anomaly even in the existence of linear divergent integral 13|, 14 ] . 
The similar verification of translational invariance can be extended to the linearly and quadratically 
divergent integrals, which is presented in the Appendix A. 

The above proof can in generally be extended to higher loops based on several theorems proved 
in ref.flU, especially based on the theorem I, theorem V and theorem VI over there. The theorem 
I is the so-called factorization theorem for overlapping divergences which states that overlapping 
divergences which contain divergences of sub-integrals and overall divergences in the general Feyn- 
man loop integrals become completely factorizable in the corresponding ILIs. The theorem V is 
the so-called reduction theorem for overlapping tensor type integrals which states that the general 
overlapping tensor type Feynman integrals of arbitrary loop graphs are eventually characterized 
by the overall one-fold tensor type ILIs of the corresponding loop graphs. This theorem is the key 
theorem for the generalization of treatments and also for the prescriptions from one loop graphs 
to arbitrary loop graphs. The theorem VI which is the so-called relation theorem for tensor and 
scalar type ILIs which states that for any fold tensor and scalar type ILIs, as long as their power 
counting dimension of the integrating loop momentum are the same, then the relations between 
the tensor and scalar type ILIs are also the same and independent of the fold number of ILIs. This 
theorem is crucial to extend the consistency conditions of gauge invariance from divergent one loop 
ILIs to higher loop ILIs. 



III. WARD IDENTITY IN MASSLESS WESS-ZUMINO MODEL 

We begin with the massless Wess-Zumino theory which is the simplest supersymmetric model. 
The Lagrangian is: 

L = -\ (d,Af - \ {d,Bf - \xH + \F 2 + \g 2 (3.1) 
F(A 2 -B 2 )+2GAB + x(A + i l5 B) X ] (3.2) 



+9 

the action is invariant, up to a total derivative, under the global supersymmetric transformation 
shown below: 

SA = ex, 5B = -ie~/ 5 x, 

5 X = -e$(A + i 75 B) + e(F + i 75 G), 

5F = e$ X , SG = -ie^x- (3.3) 

Using functional technique, one can deduce that the one-particle irreducible(lPI) Green functions 
generating functional T is invariant under the supersymmetric transformation [18]. The supersym- 
metric Ward identity we choose to check is involving two-point irreducible functions: 

<5 2 r 5 2 T 

8A(x)5A(y) 6 ~< a " ^ 6 Xa (x) Sx p (y) = ° ^ 

This could be obtained from differentiating the equantion ST = by A(x) and x(x)[19]. In the 
momentum space, we can write it as: 

T AA (p)S ja - i(^) l(3 T Xa ^ (p) = (3.5) 
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FIG. 1: Three diagrams contribute at one-loop level 



at one-loop level, Feynman diagrams contribute to this identity are shown in FIGJTJ 



r (a) _ ( P ) 



2x4g 2 



d 4 k -7^ 



(2vr) 4 k 2 
dx 



(k — p) 2 



-i8g 2 



(2^) 4 [I 2 - x(x - l)p 



dxx'j fl p ll Io(x(x -l)p 2 ) 



212 



(3.6) 



the factor 2 appears because the wave line could be A or B. And the factor 4 results from the 
fact that the fermion is a Major ana particle. We could discern this result more clearly from 
the Majorana Feynman rules given in the appendix. According to the Feynman rules we should 
calculate (xx) firstly, and then obtain the (xx) from the relation below: 



The calculation of Taa is straight forward: 



(3.7) 



r£(p) 



8g 2 



d A k 
(2tt)4 

dH 



Tr 



l^ lv {k v -p v ) 
k 2 



dx 



{k — p) 2 

2 

xp z 



(2vr) 4 
1 

dxxp 2 Iq(x(x 



Tf A (p) = 2 x 4g 2 



dH l 



[I 2 — x(x - 

i)v 2 ) - 

8 5 2 / 2 (0) 



l)p 2 } 2 

Ho) 



(3.8) 
(3.9) 



We can see immediately that the Ward identity (|3.5j) is satisfied because the integrands cancel 
out. To arrive at above results we have only carried out Dirac algebra for 7 matrices in the four 
dimensional space-time and make the shift of the integration variables. As these operations are 
all rational in a four dimensional well-defined loop regularization method, thus we conclude that 
at one-loop level the LR method indeed preserves supersymmetric Ward identity in this simple 
model. 



IV. WARD IDENTITY IN MASSIVE WESS-ZUMINO MODEL 



We are examining another supersymmetric model. The procedure is similar to what we have 
done in the above massless model. The Lagrangian of massive Wess-Zumino model is: 



L 



1 



1 



1 



1 



-F [A 2 - B 2 ) + 2GAB + x{A + %-%B)x 



(4.1) 
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It is different from the massless case with the mass term m(AF — BG — \xx)- I n this model 
bosons and fermions have equal masses as demanded by supersymmetry. In section IVIl we will 
explicitly show that the radiative corrections do not violate such an equality. The super symmetric 
transformation of component fields are the same as Eq. d3.3p . Following the same procedure, the 
two-point Ward identity of this model is extended to be[20]: 

r AA {p)s ia - i(tf) 1/3 r Xa x (p) + i(tf) 7a r AF (p) = o (4.2) 

At one-loop level, the diagrams which contribute to this supersymmetric Ward identity are 
shown in FIGJ2J FIGJ3] and FIGfJ 

It is easy to show that two diagrams in Fig.(j3]) contribute to Taf and their contributions cancel 
each other: 



?af = V J 



j£k_ 

WY 1 



1 m 1 — m 

+ 



k 2 — m 2 (k — p) 2 — m 2 k 2 — m 2 (k — p) 2 — m 



(4.3) 



The calculations of other diagrams are straightforward, we just present the final results as 
follows: 

„(a) , x , 2 [ d *k -jjfjd^ + m 1 -i^^ + m 1 

T vl?M = A 9 / 7^4 [ U2 ^2 <u ^2 ^2 +(H5) d ^2 (*75) 



Xaxew * J (2vr)4 L k 2 -m 2 (k - p) 2 - m 2 K ,0 > k 2 - m 2 v ,0 '{k-pf 



89 2 I 7^1 I dx 



d 4 l f 1 —ix^j^pV 



(2tt) 4 7 [I 2 - x(x - l)p 2 - m 2 ] 2 
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8 5 2 



dx(—ixj fl p )1 )Io(x(x — l)p 2 + m 2 



8^ 



8^ 



d 4 /fc 
(i 4 A; 



Tr 



-ij^k^ + m —i^f v (k v — p u ) + m 



k 2 — m 2 



{k — p) 2 — m 2 



(2vr) 4 Jo 



dx 



2 [m 2 + x(l — x)p 2 ] 



[I 2 — x(x — l)p 2 — m 2 ] [I 2 — x(x — l)p 2 



m 



212 



dx 



2x V 



l2(x(x — l)p + m ) — 2 m + x(l — x)p Iq{x{x — l)p + m ) (4.5) 
d 4 A: -A; 2 1 



(2^) 4 /c 2 - m 2 (A;-p) 5 



777/ 



8/ 



8^ 



d 4 / 

(2tt) 4 



dx 



-1 



+ 



777 2 + x(l — 2x)p 2 



[I 2 — x(x — l)p 2 — m 2 } [I 2 — x(x — l)p- 



77? 



212 



dx 



r<&(p) = 2 x V 



d 4 ^ 



77? 



(27r) 4 A; 2 - 77? 2 (k-pf 



dH 



dx 



m 



(2?t) 4 Jo [I 2 - x(x - l)p 2 



8g 2 / m 2 Io(x(x — l)p 2 + m 2 ) 



m 



212 



(4.4) 



-l2(x(x — l)p 2 + 777 2 ) + m 2 + x(l — 2x)p 2 Iq(x(x — l)p 2 + 777 2 ) (4.6) 



(4.7) 

Adding all the contributions together, we can see that the integrands cancel out and the super- 
symmetric Ward identity holds. Again, to arrive at above results we have only performed Dirac 
algebra for 7 matrices in the four dimensional space-time and make the shift of the integration 
variables. It further shows that in the massive Wess-Zumino model the LR method can preserve 
supersymmetry as well. 



V. WARD IDENTITY IN SUPERSYMMETRIC GAUGE THEORY 

Let us consider a more complicated case, i.e., the supersymmetric Yang-Mills theory. This model 
involves supersymmetry as well as gauge symmetry. In the Wess-Zumino gauge, the Lagrangian 
(with source terms) can be writen as: 

L = -^(F^ u ) 2 - -{d^AD + C* a dEf ab C b - -\ a IJ ab \ b + ^D 2 a + J^A^ + J a X a + f D D a (5.1) 

where A a is a Majorana spinor and D a is the auxiliary field. Similarly, the supersymmetric Ward 
identity is derived by considering the functional variation of the Green function generating func- 
tional under an infinitesimal supersymmetric transformation. All the fields transform as follows: 

5A% = -e 7M A a , 

5D a = ei^ ab X b (5.2) 



which lead to the following supersymmetric Ward identity [a. l21|: 



5 f&] 5 ^ V '\ ^V - *X- e -7^(z))7AV)> 
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(2) (3) 
FIG. 5: the 1st term 





(5) (6) 
FIG. 6: the 2nd term 





FIG. 7: the 3rd and 4th terms 



(9) 



+ 



+ 



+ 



• ^(x)^-^^ (/)> 

8A°(z) 5X b (y) 

Mfa) 5X b {y) 
5 2 J b \y' 



■{i\^y')iA c '{z'){d,Cnr f ei^ e ^)C f {x)) 



SA p (z)5X b (y) 

S 2 J b '(y') 
5AcJz)5X b (y) 



(d-A a (x)e$X a (x)iX b '(y')) 
(iX b '(y')(d p C* a )f aef H^ e (x)C f (x)) 



(5.3) 



here the notation (...) represents connected Green functions and the integration over x,y',z' are 
abbreviated. At the tree level, only the 1st, 2nd and 3rd terms in the above equation contribute, 
one can easily verify that the identity holds. At the one-loop level, only the 1st, 2nd, 3rd and 4th 
terms contribute, all the diagrams to this order are shown in FIGJ3J FIG J6] and FIG0 

We would like to point out that the first line of eq. (|5.3p is exactly the self-energy function of 
the gauge boson at one-loop level, it can be seen from the relation below: 



sA°Jz) 



r A e> Aa {z - z 

A o' A P 



(5.4) 
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Gauge symmetry requires this term to be transverse. 

We now turn to the calculation of each term in the Ward identity, and choose £ = 1 for simplicity 



n (i) 



1 



2 9 2 C acd C bdc j ' {2 ^A q 2 {q+p) 2 

-2p^p u + (5p 2 + 2p ■ q + 2q 2 )g^] 



■T,g 2 CacdCbd, 



/dx [ 
J (2tt 



d 4 l [I0y u + Wx 2 p^ Pu - 7p^p u + Ap 2 9fMU + 2{l 2 + x 2 p 2 )g, 



/ii/j 



[P - x{x - l)p 2 ] 2 



— 2^ CacdCbdc 



dx 



{lOhvv + 2<7^/ 2 ) + (10x 2 p^ " 7v»Vv + 4p 2 ^„ + 2x(2x - l)p 2 9lxv )I Q (5.5) 



n (2) 



n (3) 



9 CacdCbd, 
g 2 C ac dCbd< 



d A q 3o 



(2vr) z 



9 2 C acd C bdc J dx J 



dH 3y(P + xV) 
(2vr) 4 [* 2 - x(x - l)p 2 ] 2 



tt(1+2+3) 



35^/2 + 3x(x - l)p 2 g^h 
i r r f d^g p^qu + q^qy 

gc acd c bdc j ( ^ )4 q2(q + p)2 

f f d 4 l 1 1 

g 2 C acd C bdc J dx J (2?r)4 [f2_ a .( a ._ 1 ) p 2] 2 [-2 P ^ + V" + ^mP"'] 

g 2 C a cdCbdc J dx I 2fjLl/ + {x 2 - -^p^O 

-g 2 C acd C bcd J dx {Ax 2 - 3){p 2 g fiu - p^p v )h - AI 2fJiU + 2g flu I 2 



(5.6) 



(5.7) 
(5.8) 



Notice that in supersymmetric Yang-Mills theory, the fermions are massless and belong to the 
adjoint representation of gauge group as required by the fermion-boson symmetry. Then, 



-g 2 Atr[T a T b _ 



= g 2 C acd C bcd J dx J ^ 



d 4 q (p + q)^q u + {p + q) v q^ - {q 2 + q ■ p)g fiu 
(2tt) 4 ~" ' q 2 {q + p) 2 

{Ax 2 - 2){p 2 g^ - p^p u ) - Al^k 



[P - x{x - l)p 2 ] 2 



+ 



25, 



[I 2 - x{x - l)p 2 ] 



= g C acd C bcd / dx {Ax - 2){p 9iXV - p^p u )Io - AI 2flu + 2g /lu I 2 



(5.9) 



Adding the four terms together, we obtain the self-energy of the gauge boson which is gauge 
covariant: 



IV = g 2 C acd C bc d{p 2 g^ - p^Pu) J dxl 



(5.10) 



It is seen that the transverse condition of is satisfied in supersymmetric model with the 
Feynman gauge £ = 1 gauge. The reason is that the quadratical divergences which will potentially 
break the transverse condition cancel out in the supersymmetric model. In fact, the cancelation 
of quadratical divergences is a general feature of supersymmetric field theories, it is also one of 
the motivations to propose supersymmetry. In other word, if one wants to break supersymmetry 
but still maintain the gauge symmetry, there are several ways to realize that, for instance, give 
a mass to the fermion. In this case, the quadratical divergences do not cancel out automatically 
and they may destroy the transverse condition unless they can be regularized via an appropriate 
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regularization method to satisfy the consistency conditions [lH. As shown in 11, l^] the LR method 
is competent in this case. 

Note that here we have carried out the calculation in the Feynman gauge with £ = 1 for 
simplicity. In the general £ gauge, there is a term which could break the transverse condition if the 
regularization scheme does not satisfy the consistency condition for the logarithmic divergences, 
the term is in proportion to 

(£ - 1) * (oo - 1) 

with ao being defined via logarithmic divergent Iq^ v = ja^g^I®. In the Feynman gauge this term 
vanishes due to £ = 1. In the general £ gauge, it remains to require the regularization scheme 
satisfy the consistency condition for logarithmic divergent part, i.e., ao = 1, so that the transverse 
condition in gauge boson self-energy can hold. 

And the fermion self-energy diagram is given by: 



(2-) 4 

g 2 c acd c bcd (p 2 i a - 



CacdCbcd x 



g c acd c bcd (p 7 



(2^) 4 
dxlo 



i + 

dx 



1 



[I 2 — x(x — l)p 2 } 2 



(5.11) 



There are two diagrams from the second term of Eq. (l5.3p . The non- linear part of (FiglU(6)) 
gives rise to the contribution: 



'g 2 C ac dCbcdO' 

-g^(2p + qy] 



~n9 '' CacdCbcd 



{q+p)'- 



[g rp (p-qy + g p »(2q + p) T 



T -vh T ) 



dH 



dx 



[I 2 — x(x — l)p 



212 



:g 2 c aC dQ 



bed 



T -p z Y) / dxl 



(5.12) 



To proceed, we consider the rest diagrams coming from the third and fourth terms of Eq. (l5,3j) . 
d 4 q 



r (7) 



r<a> 



(2k) 



X [g vk (q - PY - g kfl (2q + pf + g^(2p + q) k ](-igC b dc)- 



4 [gC aC dYW + rf)7TT^H(P + g) x ^^ 



+ 



(p + qY 



-ig P k 



g 2 CacdCbcd I dx 



dH 

W) 1 



g v ^(l 2 + x 2 p 2 ) - g^p 2 - IH V - x 2 p^p u + p^p v 
[I 2 — x(x — l)p 2 ] 2 



g 2 CacdC b cd [ dx lti [g^I 2 - If + (x(2x - l)p 2 - g^p 2 + (1 - x 2 )p»p v ) J„] (5.13) 



— ^[gCdacT] u - ^ [gCbdci — f-(*g ) 



(2nT 

-g 2 CacdPbcd \ dx 



+ : 



dH 

(2«y- 



:7 



72 , ^,2^2 1„2 
2^ 



/ 2 + X V 



[Z 2 — x(x — l)p 2 ] 2 



^CacdCbcd J dxY 



h + ( x(2x - l)p 2 - -p 2 ) J 



(5.14) 
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rf) = - 9 2 c acd c bcd J 



d 4 q iq p l p 
(2tt) 



4 g 2 



+ pf 



(q + pf 



g 2 C acd C bcd I dx I — q , ^ p 



9 2 C acd C bcd J dx lp 



(2*T 

17 + 



1 V IP + x 2 p v pP - \p u pP 
[I 2 - x{x - l)p 2 } 2 



1 



p v p p h 



the total contributions of three diagrams are found to be: 

-^9 2 CacdCbcd{p 2 l u -i>p u ) J dxl 



(5.15) 



(5.16) 



After taking into account of T factors from the formula, and adding all the terms together, 
the integrands cancel out again, which demonstrates that the supersymmetric Ward identity does 
hold. To arrive this conclusion, we have only used the properties of four dimensional 7 matrices 
and translational invariance of momentum integrals. This implies that the LR method can indeed 
preserve super symmetry. The gauge symmetry holds only requiring the consistency condition for 
logarithmic divergent part due to the cancelation of quadratical divergences in supersymmetry- 
preserving regularization method. In general, to preserve gauge symmetry in non-supersymmetric 
models, it needs the consistency conditions for both quadratic and logarithmic divergences for the 
regularized ILIs. So far, we can conclude that the LR method preserves not only non-Abelian 
gauge symmetry, but also supersymmetry. 



VI. RENORMALIZATION OF MASSIVE WESS-ZUMINO MODEL 

In the previous sections we have shown that the LR method can respect supersymmetric Ward 
identities in several models including supersymmetric gauge theory, which implies that the LR 
method is viable in supersymmetric theories. While in the above applications, we have only 
used the main features of the LR method, namely the LR method is realized in four dimensions 
with translational invariance of momentum. In this section we shall apply the LR method to 
manifestly perform one-loop renormalization for the massive Wess-Zumino model. We choose 
such a model as a testing ground because it is fairly simple and well-known. The model was 
shown to be renormalizable to all orders in perturbation theory [l^| by using higher derivative 
regularization. The same conclusion can easily be obtained in the superspace formalism, where 
supergraph Feynman rules of superfields greatly simplify the calculations. For our purpose, we 
will use the component fields formalism to renormalize the theory. This is because the superspace 
formalism maintains supersymmetry in a manifest way, which is not suitable for checking the 
consistency of a specific regularization scheme in preserving supersymmetry. On the other hand, 
for the physically interesting case of broken supersymmetry, it is usually preferred to work with 
component fields. 

The action of massive Wess-Zumino model is: 

Swz = \f d*xd 2 G(^D<S> - X -m<$> 2 - X -g&) + h.c. (6.1) 

where $(x, 0, ©) is a chiral superfield. In terms of component fields the Lagrangian can be written 
as: 

L = hdpA&A + dpBdPB + ixH + F 2 + G 2 ) - m(AF + BG + ^xx) 

-g[{A 2 - B 2 )F + 2ABG + X (A - h 5 B) X ] (6.2) 
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The notions used here are slightly different from which in section IIVI It is seen that the fields 
F and G have no dynamical terms, they are auxiliary fields and can be integrated out, which is 
equivalent to eliminate them from the Lagrangian by using the equations of motions. In fact, in 
the building of phenomenological super symmetric model the auxiliary fields are eliminated by 

F = mA + g(A 2 — B 2 ) (6.3) 
G = mB + 2gAB (6.4) 

Thus the Lagrangian can be written as: 

L = ^(d^Ad^A-m 2 A 2 ) + ^(d^Bd^B-m 2 B 2 ) + ^x(^-m) X 

-mgA{A 2 + B 2 ) - gx(A - i l5 B) X - \g 2 {A 2 + B 2 ) 2 (6.5) 

which is the Lagrangian to be renormalized by using the LR method. The Lagrangian contains 
one scalar particle A, one pesudoscalar particle B and one Majorana fermion x with equal masses 
m. 

Before proceeding, we will first check what supersymmetry can tell us about the renormalization 
of massive Wess-Zumino model. The answer can easily be yielded in the superfield formalism 
based on the powerful supergraph technique. In the superfield formalism, the non-renormalization 
theorem implies that up to any order of the perturbative series only the first term (dynamical 
term) in Eq (|6.ip needs a counterterm due to the supersymmetry. Namely, after renormalization 
the action gets the following form: 

S wz = \\ d A xd 2 Q(\$>D$> - \m§ 2 - \g& + ^<5$/D$) + h.c. 
4 J 8 2 3 8 

= \ j d*xd 2 Q{\z*DS> - \^Z* 2 - \^- 2 Z^) + h.c. (6.6) 

Where the 5 term with 5 = Z — 1 is a logarithmically divergent counterterm, and Z 1 / 2 is the 
renormalization constant of the superfield. In terms of component fields, the equations of motion 
for F and G fields now become: 

F = l[ mA + g(A 2 -B 2 )} (6.7) 
G = -={mB + 2gAB) (6.8) 
After eleminating the auxiliary fields, it them leads to the renormalized Lagrangian: 
L = \z{d»Ad»A - q?A 2 ) + \z(d»Bd»B - i^) 2 B 2 ) + ±Z X (ic? - |) X 

~ ^72 ^MA 2 + B 2 ) - ^l 2 x{A - llb B) X - \(^?Z\A 2 + B 2 ) 2 (6.9) 
Which shows that the renormalizations of fields, mass and coupling constant must satisfy: 

4> bare = Z 1/2 4>; mbare = Z~ x m\ g bare = Z~ 3/2 g. (6.10) 

with (j> = A,B,x- We may summarize the features of the model: i). This model is renormalizable, 
and after renormalization all the vertexes are remained to be only one coupling constant, ii). The 
fields, mass and coupling constant share a common renormalization constant, which only contains 
logarithmical divergence. The cancellation of quadratical divergence is a general feature of all 
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supersymmetric theories, iii). As required by super symmetry, the masses of bosons still equal to 
the mass of fermion after renormalization. 

Let us now make a detailed calculation for one-loop renormalization by using the LR Method. 
The Feynman rules of Lagrangian(Eq.( |6.5p ) are listed in the appendix, there are 7 types of vertices. 
What we are going to demonstrate is that after renormalization all these 7 types of vertices will get 
the same renormalized coupling constant, and all the renormalization constants satisfy Eq. (|6.10j) . 
It is easy to verify that one- loop contributions to (A), (B), (AB), (AAB), (BBB), (AAAB), 
(ABBB) are vanishing. The rest of divergent diagrams at one-loop level are shown in FIG El the 
permutation graphs are not presented for simplicity. 




A B / V\A| fWVff B 



A WV *\AA/ y 4 



B / W\> *VWg 

.4, 





A WV *\AA/g 



FIG. 8: non-vanishing one-loop divergent graphs in massive Wess-Zumino model 

The field strength and mass renormalizations of filed A can be obtained from the calculations 
of two point Green function (AA). Five diagrams can contribute to (AA), the total contribution 
is found to be: 



L 



(AA) = \j 7^i[(- 6im 9) 2 + (-2im 5 ) 21 1 



k 2 — m 2 (k + p) 2 — m 2 



+ 



(-12ig 2 -4ig 2 



k 2 



77T 



-tr 



d A k 

^{2igC^)-^C T {2igC^)- 1 1 

(2vr) 4V y ']j,-m vy '(ty+fj-m 



C 1 
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„ 2 f 1 f d 4 k 2(1 - x)p 2 + m 2 



o 7 (2vr) 4 [fc 2 - m 2 - x(x - l)p 2 ] 2 

= A 9 2 f o J -0^12(1 -x)p 2 + m 2 ]I (6.11) 

which is only logarithmic divergent as the quadratical divergences cancel out. Using the loop 
regularization, the regularized Iq has the following explicit form: 

?' M 2 a 2 

# = T^ [ln ^-^ + yo{ wP ] (fU2) 

We shall adopt a subtraction scheme similar to the Minimal Subtraction scheme in dimensional 
regularization. For that, it is useful to introduce an arbitrary energy scale parameter /j, s and write 



l£ as: 



tf-isHf + it^ l 4-'"' +M i )] (613) 

M 2 

then the divergent terms proportional to j^^ n —§- f° r M c — ► oo in the Feynman integral are 
canceled by counterterms. As such a divergent term is independent of the Feynman parameters x, 
we can integrate x easily and obtain the divergent part of these diagrams: 

i M 2 
L{AA);div = -fTf 2 ^ 2 + rn )ln—§- (6.14) 

The counterterms correspondig to this divergence is: 



where 



from this we finally get: 



5L = U A (d^A) - \b mA m 2 A 2 (6.15) 



t 1 2 7 M* 1 2 M 2 C 

5 A = -— 2 gln-^; S mA = — 2 g ln-^ (6.16) 



where 



1 M 2 1 M 2 

Aare = (1 ~ g^^^"^)^ = ^ 1/2 A; "^A6are = (1 + 4^2#^ n ~^) m = ^ ^ ( 6 ' 17 ) 

1 2, A*? 

2 = 1 - 4^ 9 '"7# (6 ' 18) 



The calculation for (BB) is similar, which gives: 



1 iW 1 iW 

Bbare = (1 - -^g 2 ln—^)B = z 1/2 B; rriBbare = (1 + —^g 2 ln—^-)m = z~ l m (6.19) 

We now turn to the calculation of (xx)> from the Fig.® we can read directly: 

1 , /■ d 4 A; -8xg 2 Cty 
ax 



o y (2vr) 4 (A; 2 - m 2 - x(x - l)p 2 ) 2 



i M 2 

-— 2 g 2 C^ln-± + ... (6.20) 



17 



We need the following counterterms to cancel this divergence: 



where: 



1 2, M c c n 

jg In—; d mx = 



which indicates that the renormalization of field and mass is given by: 



(6.21) 



(6.22) 



Xba 



1 2 , M 2 1/2 . 1 2, M Zs -1 

(1 - ^5 ln—)x = z 1 x; m x bare = (1 + ln —) m = z m ( 6 - 23 ) 



4vr 2 ' 



So far, we have worked out the renormalization constants for the fields A, B and \ an d their 
masses, the results agree with Eq. (|6.10p . Let us switch to the renormalization of coupling constant. 
As mentioned above, there are 7 types of vertices which should be described by only one coupling 
constant when supersymmetry holds. The contributions from all divergent diagrams shown in 
FIG IS are found to be: 



L 



(AAA) 



L 



(ABB) 



L 



(A X X) 



L 



(BXX) 



-tr 
1 



d A k 

(2tt) 4 v ~ a ~ "'ty-m 

d 4 k . % 
(—6img) 



C T (2gCU 



C T (2gCU 



m +$2) - m 



C 1 



2 J (2tt) 4 



k 2 — m 2 (k + p\) 2 — m- 



■(-I2ig 2 ) + (p! -> p 2 ) + (pi -> p 3 ) 



2 7 (2vr) 4 



2img) 



k 2 — m 2 (k + p\) 2 — m 2 



■ 3 3 , M 2 . 

z — m — =- + finite terms 



-tr 



(2vr) 4 
1 /• d 4 fc 



' /1/C (2 5 Cti) ir ^-C7 T (2 5 C7t 75 ) 



\f, — m 



{-Aig 2 ) + (pi -> p 2 ) + (pi p 3 ) 



C T (2 5 Ct 75 ) 



(# -tfi)-m 



(# +i>2)-m 



(6.24) 



2 J (2vr) 4 



(—6img) 



k 2 — m 2 (k + p\) 2 — m 2 



1 

+ - 



2 J (2tt) 4 



d A k i 
{—2img) 



k 2 — m 2 (k + p\) 2 — m 2 



(-4ig 2 



(-I2ig 2 ) 



2 J (2tt) 



d k i 
\—2img) 



k 2 — m 2 (k + P2) 2 — m 2 



.1 M 2 

z — -mg in — — + finite terms 
2ir 2 



+ 



d 4 fe 

(2^F 
d 4 A; 



{2gCH 



-C T (2gCU 



(2g& 



(-4i# 2 ) + (p 2 -» p 3 ) 



C T (2 5 Ct 



(6.25) 



(27T) 4 

finite terms 

d A k , . 
(2^ 



75; 



-c T (2#ct? 



+ 



(2tt) 4 

d 4 A: 



C T (2gCh 5 ) 



m +$\)-m 



k 2 — m 
C T (2gC^ 5 ) 



k 2 



(6.26) 



(2vr) 4 
finite terms 



(2<?Ct 



+ $2) - m 



+^1) - m 



C T (2gCU 



k 2 



75; 



(# +^2) - "I 



C T (2 5 Ct 75 ) 



C T (2 5 Ct 



75, 



A: 2 — m 2 



(6.27) 
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(2vr) 4V a ; #-m v * ; (^ + ^2)-m va (# + ^2+^3)-™ 

(# +?2 +?3 +?4j - 

4 ( - 4 ^ )2 fe^ m2( ^ + pi + p2 ) 2 _ m2 + - **) + ^ - ft) 



3 Ml 
MS 

i(BBBB) = -trJ^C^'W^^^W^ 



= i—g 4 ln — + finite terms (6.28) 



7T Z 



(2vr) 4V ' 'ty-m v (#+^ 2 )-m v (^+^2+^3) -m 

C r (2gCt 75 ) J_ C T + (pa <- p 3 ) + (pa <- p 4 ) 

(# + ?2 + J»3 + #4) - m 

4(- 12 ^ 2 fc2 - m2(fc + Pl+ p 2 ) 2 - m2 + ^ 2 "> + ( p2 ^ 

4 ( - 4?g2)2 ^-^(^ + Pi+P2) 2 -m^ + fe PS) + {P2 Pi) 
3 M 2 

i—^gHn — §- + finite terms (6.29) 

d k , , % „t / n 1 



L {A ABB) = ~tr I — — i {2gCU)—^C T {2 gCU) ] C T (2gC^ 5 ) 

J (27r) 4 Ip — m (Ip+p2)-m [jp +i>2 +fv ~ m 



+$2 +^4) - m 



\JP + P2+P?,+Pa) -m 

-tr f ^^9CU)-^C T (2gCU)-—^- C T {2gCh,) 

J (27r) 4 Ip — m Vf+i>2)-m 

CT{29C] ^m + h+h+^)-m CT 

-tr [ iPh^9CU)-^C T {2gCh,) w -}- C T (2gCU) [ 

J (27r) 4 Ip-m (Ip+p3)-m (Jp +p 3 + p 2 ) ~ m 

^(2flCt 75 ) 7 -— ^— 1 C T 

W + M>2 + P3 + Pi) ~ m 



2 k 2 — m 2 (k +pi +P2) 2 — m 2 

+ h-12ig 2 )(-4ig^ ' 



2 k 2 — m 2 (k + pi + P2) 2 — fn 2 

+{ - Al92){ - Al92) \ 2 -m 2 (k + p l+ l p3) 2 -m 2 + fe -* 
1 M 2 

= i—^g 4 ln — §- + finite terms (6.30) 

(6.31) 

Introducing the following counterterms: 

5L = -dimgA^ - 5 2 mgAB 2 - 5 3 gAxx ~ hgBxi^X 

-S 5 l -g 2 A 4 -6 6 l -g 2 B 4 -5 7 g 2 A 2 B 2 (6.32) 
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with: 




(6.33) 



It is easy to check that all the renormalized vertices lead to a single renormalization constant: 



This equation, together with Eq. (l6.17p . Eq. (|6.19j) and Eq. (l6,23p . shows that the LR method works 
well in the perturbative theory of massive Wess-Zumino model. 

VII. CONCLUSION 

In this paper we have investigated the applicability of the recently developed Loop Regulariza- 
tion method in supersymmetric theories. By checking several Ward identities in various supersym- 
metric models, we have explicitly shown that the LR method is applicable to the supersymmetric 
field theories. We have also directly carried out the calculations for one-loop renormalization of 
massive Wess-Zumino model by using the LR method with string-mode regulators, the results are 
consistent with the general conclusion yielded from the supergraph technique. Once the super- 
symmetric extensions of the standard model could be discovered at the LHC, such a symmetry- 
preserving Loop Regularization method with string-mode regulators can widely be applied to the 
computations of various supersymmetric processes. 
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APPENDIX A: TRANSLATIONAL INVARIANCE OF LOOP REGULARIZATION 

The verification of translational invariance in section 2 can simply be extended to the linearly 
and quadratically divergent integrals. 

Consider firstly the quadratically divergent integral 



by rewriting the momentum factor (k — xp) 2 into (k — xp) 2 = k 2 — 2xp.k + x 2 p 2 , then replacing k 2 
by k 2 + M 2 , one has 



9bare — 



(1 + 



8vr 2 ni 



(6.34) 




(Al) 



(k - xp) 2 -> k 2 + M 2 - 2xp.k + x 2 p 2 = (k - xp) 2 + M 2 



(A2) 
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Thus the proof in the manuscript for the scalar type logarithmic loop integration can be easily 
extended to the scalar type quadratically divergent ILIs, namely 

d*k 



L 2 ^L*= lim f>f / d * (A3) 

The regularized ILIs is well-defined and allows us to shift the momentum, we then have 

L * - 1% g c " / w^kwi - ssh § c - / ft^j - * ( A4 » 

Actually, it is this translational invariance which allows us to clarify the ambiguity caused by 
the linear divergent in evaluating the triangle anomaly and CPT/Lorentz violating Chern-Simons 
term, which was shown in ref. [141 ] . To be more clear here, we demonstrate it as follows. 

Let's first present J. Jauch and F. Rohrlich's discussion on the logarithmically divergent 
integrals 23]. Considering the following integral, 

L ° - J [(jfe_p)2 +M 2]2 

and making use of the identity, 

1 1 f 1 n(a-/3)dz 



a n f3 n Jo [(a - 0)z + (5} n+l 
for ra = 2, we can rewrite the above integral as follows 



(A6) 



f d A k f 4 f 1 (p 2 -2p-k)dz _ T 

L °~J (k 2 + M 2 ) 2 2 J dk J [k 2 + M 2 + (p 2 -2p.k)zf- l0 + Lc [A7) 

The second term L c of the right-hand side is convergent, so we can safely shift the origin of k 

fyi p^z (A8) 

and the second term reads 

£ — —2 f 1 dz [ P 2 (l-^)-2p-fc 4 , . 

Lc " 2 J dZ J ik 2 + M 2 +p 2 z(l-z)]3 dk (A9) 

The term in the numerator which is odd in k will vanish. Using the identity 

k 2 ) m ~ 2 d A k _ m 2 



(k 2 + M 2 ) n (M 2 Y 



B(m,n-m) (A10) 



where B(m, n — m) = r(m)T(n — m) /T(n) and n > m > is the condition of convergence. So the 
second term in eq (IA7j) now goes as 

Lc = ~2%- fdz M f} l ; 2z) =- l7 r 2 ln[M 2 +P 2 z(l-z)]\l = (All) 
2 Jo M z +p z z(l — z) 

Therefore, for the logarithmic divergent integral, we arrive at the following identity 

f d A k f d 4 k 

L ° = J [{k-p) 2 + M 2 } 2 = J [k 2 + M 2 ] 2 = Io (A12) 
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which is independent of the regularization. 

Nevertheless, if firstly applying the Loop Regularization prescription and then shifting the 
momentum, the corresponding relation becomes a straightforward consequence 

L - Lq = lim fcf/ f. k - 1? = lim V(f/ u /\^, = (A13) 

NMff^ J [{k-p) 2 + M 2 } 2 n,m?{^ J [k 2 + M 2 ] 2 v ' 

Let us now consider the linear divergent integral. When using the identity eq. (|A6|) . a similar 
proof can be carried out and shows that a shift of A; in a linearly divergent integral will result in a 
finite additive constant 

L ^ = J l(k- P ) 2 + M 2 ] 2 = J [k 2 + M 2 \ 2 ~ — P 55 h » + + Lc ^ (AM) 

which has been shown to cause an ambiguity in evaluating the chiral anomaly if the regularization 
schemes are not applied appropriately 14]. This is because the results may depend on the procedure 
of applying the regularization schemes before or after using the identity eq. (|A6|) . 

To be safe, we shall apply LR prescription before shifting the momentum, it then leads to the 
following result 

k^d 4 k 



Ft X — ^ ./V f fed K 

Li,» - L liM = j lim J [{k _ p)2 + Mf]2 

~ NMfh J [k 2 + M 2 ] 2 ~NMfh J [k 2 + M 2 ] 2 (A15) 



where we have shifted the momentum for the well-defined regularized integral but without using 
the above identity. 

On the other hand, when applying LR prescription before shifting the momentum, but using 
the identity presented above for the integration, we then arrive at the following expression 

Wd A k 



Ll„-» Lf „ = lim Ycf [ 



[{k - p) 2 + M 2 } 2 



The second term of the right-hand side actually vanishes due to the following conditions for the 
coefficients in LR 



N 



lim Vcf(Mfr = (n = 0,l,---) (A17) 
thus we finally yield the following relation 

t r Y A N r k»d*k _ ^ N r (k+ P yd 4 k 

"™?h q J [(k-p) 2 + M 2 ] 2 NMf{^ Q Cl J [k 2 + M 2 ] 2 



N ' lvl i 1=0 



which just shows that in the LR method the translation of momentum can safely be made for a 
linearly divergent integral. 
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Now we turn to the quadratically divergent integral, 

f d A k 

L2 = J W^WTW] (A19) 

which can be rewritten as follows when using the previous identity 

f d A k f , 4 , f 1 (p 2 -2p-k)dz /A . 

L 2 = / 77o 7797 "2 d A k f-r^T, TTl? = 7 2 + L 2c A20 

J (k 2 + M 2 ) J Jo [k 2 + M 2 + {p 2 - 2p ■ k)z] 2 y ' 

Since the second term involves only linear and logarithmical divergences, we can then use the 
previous identities for those integrals when shifting the origin of k, and get the following result 
with a finite additive constant 

r - o />, f 1 {p 2 -2p-k)dz 
2c ~ J Jo [k 2 + M 2 + (p 2 - 2p ■ k)z} 2 
r 1 r P 2 d A k 
~ ~ 2 Jo J [k 2 + M 2 + P 2 z(l-z)} 2 

WW ,,., M^r^ -i^ (A21, 



[k 2 + M 2 + p 2 z(l - z)} 2 

The term which is odd in k does not contribute, and two integrals of the right-hand side cancel 
each other due to the relation 

f 1 f zp 2 d A k -I/* 1 ,/ f P 2d4k 

Jo J [k 2 + M 2 + p 2 z{l - z)] 2 ~ 2 Jo J [k 2 + M 2 + p 2 z(l - z)} 2 { ' 

Thus we arrive at the following identity 

f d A k f d A k 2 o T T /Ann . 

= J [(k- P ) 2 + M 2 ] = J wrm ~ m v L2c (A23) 

Just like the discussion in linearly divergent integral, by applying the LR prescription before shifting 
momentum, we have 

L 2 -+ L* = lim Vc^ / d J = lim V cf / * = I* (A24) 

2 n,m?,„ J (k-p) 2 + M 2 } N,M?f-t J k 2 + M 2 } 2 



i 1=0 



where the shift of momentum has been made for the regularized L^- On the other hand, again 
applying the LR prescript before shifting momentum, but using the identity obtained above, we 
arrive at the following expression 



N f d^k 



' l 1=0 



= lim Vcf / tt/^TTIT " ™ 2 P 2 lim E c ^ (A25) 

accordingly, because of the vanish of the second term in the right-hand side, we obtain the same 
regularized result 

N r rl 4 k N r d A k 

L§ = lim Vcf / Wl -g. = lim VcW n9 aj2l = I? (A26) 

N > M ?T=1) J [{k-p) 2 + M 2 } N,M?f^ Q 1 J [k 2 + M 2 ] 2 v ; 

So far we have demonstrated that Loop Regularization can preserve translational invariance not 
only in logarithmically, but also in linearly and quadratically divergent integral. 
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APPENDIX B: DERIVATION OF MAJORANA FEYNMAN RULES 



Here we are going to present a simple and definite derivation of Majorana Feynman rules which 
are useful for our calculations in this paper. We will begin with the quantization of free Majorana 
fermion, and figure out the difficulties of formulating the Majorana Feynman rules, then provide a 
consistent prescription. The unusual Majorana Feynman rules are result from the Majorana fermion 
self-conjugacy. Though the two-components formulation of Majorana field is more fundamental, it 
is still very useful to work in four-components formalism because the 7 matrices is more convenient 
for practical calculations. 

The Majorana fermion field x is quantized by stipulating the following equal-time anticommu- 
tators: 

{Xa(x),Xp(y)} = <W<5 3 (x-y) 

{xa(x), X/3 (y)} = {xl(x),4(y)} = o (Bl) 

The plane wave decomposition of x is n °t obvious. In two-components formalism the difficulty 
behaves as that the equation of motion (EOM) is no longer a linear equation since the EOM 
connects \ to its complex conjugation. In four-components formalism the difficulty lies in the 
Majorana condition: x = X° = Cx T ■ But if we use the spinors u and v which satisfy s = Cv^ g 

and ujj. g = Cu^ s , then x can be expanded as: 

x = I Wf i ^ [c ^V lkx + c isV ikx ] (B2) 

here c and are the annihilation and creation operators of Majorana fermions. For Majorana 
fields, we still have: 

(0\Txa(x)Uy)\0) = J ^ e ~ lk(x ~ y) (j^U = - y) (B3) 

Note that because of the Majorana condition x = Cx T and % = X T C, (0\Txa(%)xp(y)\fy an d 
{0\Tx a (x)xf3(y)\0) do not vanish. It is easy to show that: 

(0\T Xa (x)xp(y)\0) = S Fai {x-y)C^ (B4) 
(0\Tx a (x)xf3(y)\0) = C^S Fj0 (x-y) (B5) 

The explicit expressions of u-^ and s as well as the spin-sum identities can be found in j22|, we 
list the results here: 

«k„ = (^;). Sk ,.=(cjv^, c «vw) 

= ( S ) ' ^ = (- 2 <J^ 2 <JV^) <*» 
and C±i/2 are defined as below (here 6 is the polar angle of k, and 4> is the azimuthal angle of k.): 
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The spin-sum identities are: 



\ 111 

s 


C4 \r 

K,s 


— // —1— 771 


\ lit 

s 




— — 771 

/X III/ 


s 


T 


= [fjt + m)C 


I>k,, 




= m-m)C T 




^k, s 


= C\lf-m) 




^k, s 


= C\lf + m) 



(B8) 



Before starting the derivation of Majorana Feynman rules we may briefly review the derivation 



for the usual Dirac fermions. The argument below follows the one in 2J]. The calculation of a 



typical scattering matrix element corresponds to the evaluation of the following expression: 

(0|6 1 ...6 m d 1 ..AT[(^(x 1 )r^(x 1 ))...(^(x / )r^(xO)]6l..^4...4|0) (B9) 

Firstly, we should rearrange the interaction terms to make them following the order of contractions. 
Since only one type of contraction (ipijj) exists for Dirac fermion, the internal propagator reads: 
(t/jif)) = Sf{p), here the fermion charge and the momentum flows are well defined from ip to ip, 
the Feynman rule for vertex directly reads as ir. For Dirac fermion, the fermion charge flow (in 
fact this is also the momentum flow) of internal popagator forms a continuous flow, when writing 
down the analytic expression one should first do it oppositing to the continuous flow. The most 
important step is to determine the Relative Sign of Interfering Feynman graphs (RSIF). There are 
in general three types of commutations which can contribute to the RSIF. Firstly, when reordering 
bi, di, b\ and d\ to put them in the appropriate places of Wick contractions, it causes a factor 
(— l) p . Here P is the parity of the permutation of the annihilation and creation operators. This 
factor can be read from the order of external spinors in the analytic expression with respect to 
the given reference order. Secondly, for a closed fermion loop, one needs to exchange the first and 
the last field operator in the fermion chain, which gives a factor (— 1) L , where L is the number 
of fermion loops. Finally, since d\ must contract with ip and di must contract with ip, one needs 
to move the creation operator di to the beginning of Wick contraction and move the annihilation 
operator d\ to the end, which leads to a factor (—1)^ with V being the total number of spinors v 
and v. Since V is universal for all graphs of a given process, this factor can therefore be omitted. 

We now trun to investigate the Majorana fermion case. Firstly, we consider the situation 
that there are no Dirac fermions but only Majorana fermions. As mentioned above, all possible 
contractions between x and/or \ do not vanish now. In this case, after rearranging the interaction 
terms to perform Wick contraction for operators one by one, we need to consider four types of 
Majorana propagators, i.e. (xx)> {xx)i (xx) an d {xx}- More seriously, the propagators depend on 
the sign of its momentum p, but now we can not define the orientation from x to % as the arrow 
of momentum. That means we need to find out a new method to resign the arrow of momentum. 
For the Feynman rule of vertex, it raises a new ambiguity. For instance, when contracting an 
interaction Lagrangian xTx m the- time-order product, one can contract the operator x with one 
field operator lies on the left of this vertex and contract x with another lies on the right, or one 
can also contract x with one field operator lies on the left and contract x with another lies on 
the right. In the later case an additional (—1) will emerge. Previous discussions [25j, [26| for the 
Majorana Feynman rules follows this analysis and try to reduce the number of propagators and 



25 



vertices, while the resulting consequences are still too obscure and not easy to use. In ref 24], 
the author introduced the charge-conjugate fields ip c and ip c to Feynman rules and tried to give a 
uniform description of Dirac and Majorana field. Here we shall provide an alternative and simple 
description. 

Firstly, we may eliminate x from the interaction Lagrangian by using the Majorana condition 
X = ~x T C\ so that only one type of propagator (xx) remains. We then use a line without arrow 
to represent a Majorana propagator. Since Majorana fermions can not carry any charge, this 
representation is natural. In the momentum space, the Feynman rule for Majorana propagator is 
j^^C T . To obtain the Feynman rule of vertex, we may rewrite xTx as: 

XcFapXp = Xa(-Cl p T pf3 )xf3 = -X^(rJ /0 Cj a )x Q = - Xa (-C^T - T T C f ) a(3 xi3 

lx- T 'apXp (BIO) 



2 



with: 



r' = -c f r - r T c t = -r ,T (bh) 



Now the ambiguity mentioned about disappears as V is antisymmetric. The Feynman rule for 
vertex simply becomes: iV. One can treat the Majorana fermions just like a real scalar boson to 
obtain the correct symmetric factor of a given graph. 

Next, we should determine the direction of momentum in Majorana propagators. Remember 
that generally a factor e~ lkx means momentum k flows in the point x and e lkx means momentum 
k flows out the point x. Every contraction between two field operators 0(x), 0{y) can always be 
written in the form: (O(x)O(y)) = J e -iK x ~n) S (k) , for example in our case: 

(o\T Xa (x) X p(y)\o) = J ^ e ' lk(x ~ y) (j^i cT U (B12) 

which indicates that the direction of momentum flow is always opposite to the direction of con- 
traction for a propagator, and in a fermion chain the momentum flows of propagators form a 



continuous flow its direction is opposite to the direction of contractions. In [24j] such a folw was 
called as 'fermion flow', here we may, more precisely, call it as 'fermion momentum flow'. This 
comes to the conclusion: for each fermion chain we fix an arbitrary orientation (fermion momentum 
flow), the momentums of all fermion propagators follow this orientation, and we should write down 
the Feynman rules proceeding opposite to the chosen orientation. 

Finally, to complete the Majorana Feynman rules, it needs to give the rules of external fermion 
lines and determine the RSIF. The rules of external fermion lines can easily be obtained from the 
plan- wave decomposition of %, see Eqn. (|B2p . Since 



<0|c kiS X a (x) v ak ^ k * (B13) 
Xa (x)cljO) — u^e-** (B14) 

which implies that the creation of a Majorana fermion corresponds to a spinor u ^ with momen- 
tum k flow out, and the annihilation of a Majorana fermion corresponds to a spinor u ^ g with 
momentum k flow in. If the spinor locates at the beginning of contraction, we should write it as a 
row vector say a u T or v T . Now we can give a prescription to fix RSIF. Factor (— l) p can be got 
from the permutation parity of the spinors in the obtained analytical expression with respect to 
some reference order. Factor (— 1) L can be got from the number of closed fermion loops. Factor 
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(— 1) now is a little different from which in Dirac field theory. Since moving any one creation op- 
erator arising from the initial state to the beginning of contraction will contribute a factor —1, and 
moving any one annihilation operator arising from the final state to the end of contraction which 
also contribute a factor —1, it seems that we should count the total number of such operation. 
Suppose that there are 'a' fermions in the initial sate and '6' fermions in the final state, and we 
must move ith fermion creation operators to the beginning and jth fermion annihilation operators 
to the end, then we have a — i + j = b + i — j, i.e. \i — j\ = ||a — b\. Namely, V = ||a — 6|. Since 
a and b is universal for all graphs of a process, we can always omit (—l) v all the time. 

Let us consider the situation that a Majorana fermion \ couples to a Dirac fermion ip, the 
interaction Lagrangian contains the following terms: 

xTip + i)tx {where : f = 7°r t 7°) (B15) 

When keeping a continuous "fermion momentum flow" for a fermion internal line, we then need 
to consider two types of Dirac propagators: (tpip) and {tpip) which have the following explicit forms: 

(0|T^(x)^(y)|0) = J ^ e -^-y) ( _±_ U (m6 ) 

<0|r^ a (x)^(y)|0) = J ^ e ~ lk{X ' y) i(j^) T U (B17) 

We then need to use a line with arrow to represent the Dirac propagator, the arrow reflects the 
flow of charge which flows out of tp and into tp. If the direction of charge flow coincide with the 
direction of the "fermion momentum flow", we should use {ipip) = jz^, otherwise we should use 

The Feynman rules for vertexes are also doubled. For the Dirac-Dirac interaction, one has: 

-tpaFapIpp = r 4>a{-T T ) a p'4)p (B18) 

If the direction of charge flow coincide with the direction of the momentum flow, we should use 
iT, otherwise we should use — iT T . the vertexes rules of Majorana-Dirac interaction can be derived 
similarly from the identities: 

XaFapipp = Xa{-tfT) a pil)p = ip a (-r T Ci) a pXp (B19) 

Tpa^apXp = Xa{-F T )af3i>f3 (B20) 

The RSIF can be determined by using the same method as we mentioned above. 

With the above considerations, we can summarize our Feynman rules. The solid lines are still 
used to denote the fermions. Dirac fermions lines carry arrows which reflect the direction of charge 
flow, Majorana lines do not carry arrows. We may write down Feynman amplitudes according to 
the following steps: 

1. Draw all topologically distinctive, connect Feynman diagrams for a given process. 

2. Fix an arbitrary direction for each fermion chain. This is the direction of "fermion momentum 
flow" , which means that the momentum of every internal fermion line should follow this direction. 
We should write down the Dirac matrices proceeding opposite to the chosen direction through the 
chain. 

3. For the external fermion lines, the rules are shown in FIGJ9J 

If the spinors locates at the beginning(end) of contraction, we should add a superscript T 
appropriately to write them as row(column) vectors. 

4. For the fermion propagators, the rules are shown in FIGI101 
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Creation of a Majorana fermion : * V Us Annihilation of a Majorana fermion : ' * U k s 

* P P 

Creation of a Dime fermion : " M k,s Annihilation of a Dirac fermion : '. ~ * W k,s 

•" p J J ► p 

Creation of a Dirac anti — fermion : * * ^k.s Annihilation of a Dirac anti — fermion : " * k.s 

« p ► p 

FIG. 9: Feynman rules for external lines 



Majorana propagator : fi — m 

P 

. 1 . (—^--7 

Dirac propagator : p 1 — m p" + rn 

P P 

FIG. 10: Feynman rules for propagators 



5. For the general fermion interactions xTix, $T2ip, + ^P^3X^ where T3 = 7°r37°, the 
Feynman rules are shown in FIG 1 1 1 1 respectively. 

6. To determine the RSIF. For each diagram, multipling by a factor (—1) for each closed fermion 
loop, and multipling by the permutation parity of the spinors in the obtained analytical expression 
with respect to some reference order. 

7. Multipling a symmetry factor S^ 1 for each diagram. The Majorana fermions may be treated 
just as real scalar fields to obtain the symmetry factor. 

S = g [] 2 /3 (n!) Q " (B21) 

n=2,3,... 

where a n is the number of pairs of vertices connected by n identical self-conjugate lines, is the 
number of lines connecting a vertex with itself, and g is the number of permutations of vertices 
which leave the diagram unchanged with fixed external lines. 
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i> x ir 3 x 1 V> -«r 3 T 
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FIG. 11: Feynman rules for vertexes 

APPENDIX C: FEYNMAN RULES OF MASSIVE WESS-ZUMINO MODEL 

We present all the Feynman rules of this model in FIGI121 
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FIG. 12: Feynman rules of massive Wess-Zumino model 
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